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Abstract 

We study renormalized quenched strong-couphng QED in four dimensions 
in arbitrary covariant gauge. Above the critical coupling leading to dynam- 
ical chiral symmetry breaking, we show that there is no finite chiral limit. 
This behaviour is found to be independent of the detailed choice of photon- 
fermion proper vertex in the Dyson-Schwinger equation formalism, provided 
that the vertex is consistent with the Ward-Takahashi identity and multiplica- 
tive renormalizability. We show that the finite solutions previously reported 
lie in an unphysical regime of the theory with multiple solutions and ultravi- 
olet oscillations in the mass functions. This study supports the assertion that 
in four dimensions strong coupling QED does not have a continuum limit in 
the conventional sense. 



I. INTRODUCTION 

A useful approach to studying the mechanism of dynamical chiral symmetry breaking 
(DCSB) is through the Dyson-Schwinger equation (DSE) formalism Ijl]-^. The infinite set 
of coupled DSE's must always be truncated at some point, but we can still make progress 
by closing off the tower of equations with a suitable Ansatz consistent with all appropriate 
symmetries of the theory and having the correct perturbative limit. While not a complete 
first principles treatment of a theory, this approach is nonetheless a useful tool, since it 
does allow Lorentz covariance to be maintained as well as allowing for the infrared (IR) 
and ultraviolet (UV) limits to be numerically taken in a straightforward way. In addition, 
model independent results following from symmetry priciples alone can still be obtained and 
numerically verified in a rigorous way as we will see. On the other hand, lattice gauge theory 
(LGT) studies Q are a first principles, approximation-free approach, but they still present a 
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significant computational challenge when attempting to verify ultraviolet and infrared limits 

m- 

The abelian nature of quantum electrodynamics (QED) in many ways makes it a much 
simpler system to study than a nonabelian theory such as quantum chromodynamics (QCD). 
For this reason it has been the subject of many nonperturbative studies, which have as their 
long-term goal a detailed understanding of nonperturbative QCD. On the other hand, strong- 
coupling QED4 is widely anticipated to behave unconventionally in the continuum limit and 
for this reason is a theory of considerable interest in its own right. 

In previous work |^ we first introduced a numerical renormalization procedure and ap- 
plied it to QED4 with a quenched photon propagator using the Dyson-Schwinger formalism. 
This is a direct application of the standard renormalization procedure to the nonlinear self- 
consistent framework needed to study dynamical chiral symmetry breaking. This initial 
work, in Landau gauge, was recently generalized to arbitrary covariant gauges 0. The 
central result of these two works was to demonstrate that the numerical renormalization 
procedure works extremely well and allows the continuum limit (A 00) to be taken nu- 
merically, while giving rise to stable finite solutions for the renormalized fermion propagator. 

In this article we investigate the chiral limit in renormalized quenched strong-coupling 
QED4, using a photon- fermion vertex that satisfies the Ward-Takahashi Identity (WTI) and 
makes the fermion DSE multiplicatively renormalizable. We find that for couplings above the 
chiral critical coupling, keeping the bare mass mo (A) = as the cutoff is relaxed results in a 
dynamical mass function with no oscillations and which diverges with the cutoff. The finite 
solutions described in previous articles ||^,|] showed damped oscillations in the dynamical 
mass functions at large p^, which suggested that they were unphysical. Further, we show 
here that for a given supercritical coupling and the same bare mass mo(A), it is possible to 
have multiple solutions corresponding to different renormalized masses m{fi). We conclude 
that quenched strong-coupling QED in four dimensions does not have a chiral limit in the 
conventional sense above the chiral phase transition. 

In Sec. m we briefly summarize the renormalized Dyson-Schwinger equation formalism, 
the numerical renormalization procedure, and a particular fermion-photon vertex ansatz 
that we use to illustrate our general arguments numerically. In Sec. |T| we demonstrate the 
scaling of the DSE solutions with zero bare mass for supercritical coupling, and the existence 
of multiple solutions with the same bare mass; we also present a general argument to show 
that a vertex which is consistent with the WTI and multiplicative renormalizability leads to 
a diverging mass function above critical coupling in the continuum limit. We discuss these 
results and their relevance to QCD and to unquenched QED4, in Sec. |V[ For further details 
and references and an expanded discussion we refer the reader to Refs. 



II. FORMALISM 

Dynamical chiral symmetry breaking (DCSB) occurs when the fermion propagator de- 
velops a nonzero scalar self-energy in the absence of an explicit chiral symmetry breaking 
(ECSB) fermion mass. We will refer to coupling constants strong enough to induce DCSB 
as supercritical and those weaker as subcritical. We write the fermion propagator as 
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Siv) = ^^^'^ = 1 fl) 

where we refer to Alp"^) = 1/Z{p'^) as the finite momentum-dependent fermion renormahza- 
tion and where M{p'^) = BijP') / A{p^) is the fermion mass function. In the massless theory 
(i.e., in the absence of an ECSB bare fermion mass mo(A)) by definition DCSB occurs when 
M(p2) ^ 0. 

With the exception of Refs. 0,0, most studies have neglected the issue of the subtrac- 
tive renormahzation of the DSE for the fermion propagator. Typically these studies have 
assumed an initially massless theory and have renormalized at the ultraviolet cutoff of the 
loop integration, taking Zi = Z2 = 1- Where a nonzero bare mass has been used, it has 
simply been added to the scalar term in the propagator. Although there were earlier formal 
discussions of renormahzation p|,|7|-|l0[] , the subtractive renormahzation program had not 
previously been implemented. 

We will concentrate our discussion on quenched strong-coupling QED4, where here the 
term "quenched" means that the bare photon propagator is used in the fermion self-energy 
DSE, so that Z3 = 1 and there is no renormahzation of the electron charge. It should be 
carefully noted that this is a slightly different usage from that found in lattice gauge theory 
studies, since in DSE studies with a quenched photon propagator virtual fermion loops may 
still be present in the proper fermion-photon vertex. 

The DSE for the renormalized fermion propagator, in an arbitrary covariant gauge, is 

S-\p) = Z,{f^,A)[^~mo{A)]-tZ,{f^,A)e' f -0^^YS{k)r%k,p)D,,{q) ; (2) 

here q = k — p is the photon momentum, /i is the renormahzation point, and A is a regular- 
izing parameter (taken here to be an ultraviolet momentum cutoff). We write mo(A) for the 
regularization-parameter dependent bare mass. The renormalized charge is e (as opposed 
to the bare charge Cq), and the general form for the renormalized photon propagator is 

with ^ the renormalized covariant gauge parameter and ^0 = Z3{fi,A)C, the corresponding 
bare one. In the quenched approximation, we have for the coupling strength and gauge 
parameter respectively a = e"^ /in = ao = el/in and ^ = ^o; and for the photon propagator 
we have 

Z^-(g) ^ Di;^iq) = I + ^) - ^ • (4) 

The requirement of gauge invariance in QED leads to a set of identities referred to as 
the Ward-Takahashi Identities (WTI). The WTI for the fermion-photon vertex is 

q,T^{k,p) = S-\k)-S-\p) , (5) 

where q = k — p . This is a generalization of the original, differential Ward identity, which 
expresses the effect of inserting a zero-momentum photon vertex into the fermion propagator. 
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dS-^{p) 
dpu 



r^(p,p). 



(6) 



The Ward identity, Eq. (^), follows immediately from the WTI of Eq. (|^). In general, for 
nonvanishing photon momentum q, only the longitudinal component of the proper vertex 
is constrained, i.e., the WTI provides no information on ri^{k,p) = T^^T^ip^k) for g 7^ 0. 
[We use the notation T^'^ = g^" — {q^q^ /q^) and C^'^ = {q^q'^ /q^) for the transverse and 
longitudinal projectors, respectively.] In particular, the WTI guarantees the equality of 
the propagator and vertex renormalization constants, Z2 = Zi (at least in any sensible 
choice of subtraction scheme The WTI can be shown to be satisfied order-by-order in 
perturbation theory and can also be derived nonperturbatively. 

As discussed in []T],|rT|, this can be thought of as just one of a set of six general requirements 
on the vertex: (i) the vertex must satisfy the WTI; (ii) it should contain no kinematic 
singularities; (iii) it should transform under charge conjugation (C), parity inversion (P), 
and time reversal (T) in the same way as the bare vertex, e.g.. 



C-'T,ik,p)C = -TU-p,-k) 



(7) 



(where the superscript T indicates the transpose); (iv) it should reduce to the bare vertex in 
the weak-coupling limit; (v) it should ensure multiplicative renormalizability of the DSE in 
Eq. (^; (vi) the transverse part of the vertex should be specified to ensure gauge-covariance 
of the DSE. 

Ball and Chiu |T^ have given a description of the most general fermion-photon vertex 
that satisfies the WTI; it consists of a longitudinally-constrained (i.e., "Ball-Chiu") part 
EgQ, which is a minimal solution of the WTI with no artificial kinematic singularities, 
and a basis set of eight transverse vectors Tf(A;,p), which span the hyperplane specified 
by Cfj_yT^{k,p) = (i.e., qyT^{k,p) = 0), where q = k — p. The minimal longitudinally 
constrained part of the vertex will be referred to as the Ball-Chiu vertex and is given by 



^Bcik,p) 



A{p')]r + 



{k + pY 

^2 _ p2 



[Aik')-Aip% 



[B{P) - Bip')] 



(8) 



Note that since neither C^u^'^ci^^P) '^fj.i^^Bci^^P) vanish identically, the Ball-Chiu ver- 
tex has both longitudinal and transverse components. The transverse vectors can be conve- 
niently written as [|l^ 



Tr{k,p 
T^{k,p 

n{k,p 

mk,p 
T^k.P 
T^k^P 

mk,p 

t^Kp 



p^{k ■ q) - k^{p ■ q) , 
[p^'ik ■ q) - k^ip ■ q)] 



^-^-q^4.. 



q J 
,2 



Y{p^ - k^ 
1 

2 

i^k'^p^a^x + k^'^ 



A;^] - 2i{p - kYk^p^axv 



+ {p + kY4, 
(p^ -k^)[Yi^+ j^)-p'' 



k'']-t{k + p)''kY(^x. 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 
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where we use the conventions g^'^ = diag(l, — 1, — 1, — 1), {7'*,7'^} = ^9'^'^, a-nd cr'^'^ = 
(i/2)[7'*, 7^]. A general vertex is then written as 

r'^{k,p) = T%c{k,p) + i:n{k^p^ q')Tt{k,p) , (17) 

1=1 

where the are functions that must be chosen to give the correct C, P, and T invariance 
properties. 

As previously mentioned, the renormalization procedure is entirely standard. One first 
determines a finite, regularized self-energy, which depends on both a rcgularization parameter 
and the renormalization point. One then performs a subtraction at the renormalization 
point, in order to define the renormalization parameters Zi, Z2, Z3 which give the full 
(renormahzed) theory in terms of the regularized calculation. Consider the regularized self- 
energy E'(//, A;p), leading to the DSE for the renormalized fermion propagator, 

S-\p) = Z2(/x, A)[^- mo(A)] - A;p) 

= ^-mii^) - = A{p^) fi-B{p^) , (18) 

where Tj{ix;p) denotes the renormalized self-energy and the regularized self-energy is given 
by (g = (k-p)) 

E'(/., A;p) = zZi(/., A)e' f -^,l^S{^x■ k)T^{^x■ k,p)D^,{i,; q) . (19) 

Here D^'^{ij;q) and r'^{iJ,]k,p) denote the renormalized photon propagator and photon- 
fermion proper vertex respectively. As suggested by the notation (i.e., the omission of 
the A-dependence) renormalized quantities must become independent of the regularization- 
parameter as the regularization is removed (i.e., as A — > 00) in a renormalizable theory. The 
self-energies are decomposed into Dirac and scalar parts, 

E'(//, A;p) = E'.ii,, A;p') i> + E;(//, A;p^) (20) 

(and similarly for the renormalized quantity, E(/x,p)). By imposing the renormalization 
boundary condition, 

S-\p) ,=^-m(/x), (21) 

one gets the relations 

E,,,(/i;p2) = E:,_,(//, A;p^) - E:,,,(//, A; /x^) (22) 

for the self-energy, 

Z2(//,A) = 1 + E:,(//,A;//2) (23) 
for the renormalization constant, and 

mo(A) = \m{ix) - T.'X^x, A; /x^)] /^^(/x. A) (24) 
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for the bare mass. The mass renormahzation constant is then defined as 



A) = mo(A)/m(/i) 



(25) 



i.e., as the ratio of the bare to renormahzed mass. The vertex renormahzation, A), is 

identical to Z2(/i, A) as long as the vertex Ansatz satisfies the Ward Identity; this is how it 
is recovered for multiphcation into E'(/i, A;p) in Eq. (|19|). 

The chiral limit occurs by definition when the bare mass is set to zero and the regulariza- 
tion is removed, i.e., maintaining mo (A) = while taking the limit A oo. Explicit chiral 
symmetry breaking (ECSB) occurs when the bare mass mo (A) is not zero. Dynamical mass 
generation or dynamical chiral symmetry breaking (DCSB) is said to have occurred when 
M{p^) 7^ in the absence of ECSB. As the coupling strength increases from zero there is 
a transition to a DCSB phase at the critical coupling strength a^- Concisely, the absence 
of ECSB means that mo (A) = and the absence of both ECSB and DCSB (i.e., a < Oc) 
means that M{p'^), '"^(yu), and mo (A) simultaneously vanish. (Recall that in the notation 
that we use here, M{p'^) = B{p'^)/A{p'^) and m(/i) = M(/i^).) This is the same definition 
of the chiral limit that is used in nonperturbative studies of QCD, see e.g. Refs. and 
references therein. Obviously, any limiting procedure where we take mo (A) sufficiently 
rapidly as A ^ oo will also lead to the chiral limit. 

Let us temporarily indicate explicitly the choice of renormahzation point by a /i- 
dependence of the renormahzed quantities, i.e., A^fi-jp"^) = l/Z(/i;p^), M(/i;p^) = 
B{fi;p'^)/A{^;p'^), etc. Note that Eq. (IIS) implies that 



A{i^;p') 
B{^,;p') 



Z2(/i, A) - S:,(/i, A; = 1 - Srf(/i, A; p^) , 
Z2(/i,A)mo(A) + S;(/i,A;p2) = m(/i) + S,(/i, A; 



(26) 



The renormahzation point boundary condition in Eq. ( ^I]) then leads to S(/i, A; /i^) = 0, or 
equivalently, to the two boundary conditions A(/i; /i^) = 1 and M(/i; /x^) = -B(/i; /i^) = m(/i). 
From Eq. (E^) we have 



A(/i;p2)/Z2(/i,A) 



1 



S:,(/i,A;p2)/Z2(^,A) 



Bifi;p'')/Z2ifx,A) =mo(A)+ A;p^) /Z^ifi, A) 



(27) 



The renormahzation group is the set of renormahzation point transformations which by def- 
inition leave the bare quantities of the theory unchanged. Hence, since mo (A) is renormal- 
ization point independent it is clear from Eq. ( jn\ ) that y4(/i;p^)/Z2(/i, A), i?(/i;p^)/Z2(/i, A), 
S^(/i, A;p^)/Z2(/U, A), and E^(/i, A;p^)/Z2(/i, A) are renormahzation point independent. 
Hence S{n;p)Z2{n,A) is renormahzation point independent. In shorthand form we can ex- 
press this as S{fi;p) oc l/Z2(/i, A) under a renormahzation point transformation. Hence, the 
choice of renormahzation point is equivalent to the choice of scale for the functions A and B. 
Similarly, in the general unquenched case |T^ under a renormahzation point transformation 

we have in addition D°''^(n; q) oc ^(/x) oc l/Z3(/i, A), e(/i) oc 2'2(/i, A) y^Za (/x. A) /Zi (/i. A), and 
Vdu,] q,p) oc Zi(/i, A). It is straightforward to verify the consistency of these renormahza- 
tion point transformations. For example, from Eq. (|19|) we see that these scaling properties 
ensure that S' oc Z2 and hence from Eq. ( pTsD that S{p) oc I/Z2 as it should. Thus, since 
Zi = Z2, multiplicative renormalizability will automatically follow if we ensure that F — cF 
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as A{p'^) cA{p'^) and B{p^) cB{p'^). This behavior is automatic for the Ball-Chiu part 
of the proper vertex rgQ(/i;p, A;) as can be seen from Eq. (§). This consistency is unaf- 
fected by considering the quenched photon propagator case where = 1. Clearly in order 
to ensure multiplicative renormalizability for an arbitrary vertex it is necessary to choose 
the functions to scale in the same way, i.e., r, oc CTj. This is the precise statement of 
the restriction that multiplicative renormalizability imposes on the proper photon-fermion 
vertex. 

From the above arguments we see that under a renormalization point transformation we 
must have for all 



M{fi';p^) = M{fi;p^) = M{p') , 

A{^^;p^) Z2(/i,A) ""^^'^^ ^'^^ 

from which it follows for the fermion propagator that S{fi']p)/ S{fi]p) = Z2(/i, A)/Z2(/i', A) 
in the usual way. The behavior in Eq. ( ]28|) is explicitly tested for our numerical solutions. 
It is clear from Eq. (^Sf) that having a solution at one renormalization point (/i) completely 
determines the solution at any other renormalization point {fi') without the need for any 
further calculation. 



Example vertex choice 

Our results and conclusions regarding the chiral limit are general and do not depend on 
any specific vertex choice, i.e., any specific choice for the functions Xj. We only require that 
the vertex satisfy the Ward-Takahashi identity and that it be consistent with multipicative 
renormalizability. However, in order to discuss the renormalized finite solutions it is neces- 
sary to present detailed numerical calculations. For this purpose we will use as an example 
the modified treatment of the Curtis- Pennington vertex introduced in Ref. 

Curtis and Pennington published a series of articles 0"TD] describing their specification 



of a particular transverse vertex term, in an attempt to produce gauge-covariant and mul- 
tiplicatively renormalizable solutions to the DSE. In the framework of massless QED4, they 
eliminated the four transverse vectors which are Dirac-even and must generate a scalar term. 
By requiring that the vertex r^{k,p) reduce to the leading log result for k ^ p they were 
led to eliminate all the transverse basis vectors except Tq, with a dynamic coefficient chosen 
to make the DSE multiplicatively renormalizable. This coefficient had the form 

Mk',p',q') = -^[Aik')-Aip')ydik,p) , (29) 

where d{k,p) is a symmetric, singularity-free function of k and p, with the limiting behavior 
\imk2 d{k,p) = k"^. [Here, A(p^) = l/Z^p"^) is their l/jF(j9^).] For purely massless QED, 
they found a suitable form, d{k,p) = {k"^ — p^Y /{kP' This was generalized to the case 

with a dynamical mass M{p^), to give 

,(,^,^j]^^j^nm^±^£^, (30) 
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It is clear that this choice scales as described in our discussion above and is consistent with 
multiplicative renormalizability. 

One refinement in our application of the C-P vertex |^ in the present work is associated 
with subtleties in the ultraviolet regularization scheme. Although there have been some 



exploratory studies of dimensional regularization for the DSE |T^, this has not yet proven 
practical in nonperturbative field theory and momentum cutoffs for now remain the regular- 
ization scheme of choice in such studies. Naive imposition of a momentum cutoff destroys the 
gauge covariance of the DSE because the fermion self-energy integral contains terms, related 
to the vertex WTI, which should vanish but which are nonzero when integrated under cutoff 



regularization [|T^,|T^]. Based on these considerations a "gauge-covariance-improved" treat- 
ment of the C-P vertex was proposed, which consists of the replacement (in the quenched 
approximation) S'^(/x, A;p^) -|- A)a^/87r — >• S^(/i, A; j9^). Full details and the derivation 
of this modification can be found in Appendix A of Ref . . This quenched approximation 
correction does not spoil the scaling needed for multiplicative renormalization since in the 
quenched approximation the correction term scales with Z2(/i, A) as does E'. We have now 
fully specified the sample vertex that we use in our numerical calculations. 

For the numerical calculations the equations are separated into a Dirac-odd part 

describing the finite propagator renormalization A(p^), and a Dirac-even part for the scalar 
self-energy, by taking ^Tr of the DSE multiplied by ^/p^ and 1, respectively. The equations 
are solved in Euclidean space and so the volume integrals, / d'^k, can be separated into angle 
integrals and an integral / dk"^; the angle integrals are easy to perform analytically, yielding 
the two equations which will be solved numerically. 

In order to obtain numerical solutions, the final Minkowski-space integral equations are 
first rotated to Euclidean space. They are then solved by iteration on a logarithmic grid 
from an initial guess. The solutions are confirmed to be independent of the initial guess 
and are solved with a wide range of cutoffs (A), renormalization points (/i), couplings [a), 
covariant gauge choices (^(/x)), and renormalized masses {m{fi)). As has been reported in 
detail elsewhere choosing the renormalized mass, m{fi), and then solving for the bare 
mass, mo(A), leads to extremely well-behaved finite solutions for v4(p^) and M(p^), which 
do not vary as we take the continuum limit (A oo). In addition, the renormalization 
point transformation properties given in Eq. (|28|) have also been explicitly verified using our 
numerical solutions, which typically have an accuracy of better than 1 in 10^. 

An important outcome of these studies was the observation that all solutions which were 
well-behaved in the continuum limit contained decaying oscillations in the mass function 
MijP') if one looked sufficiently far into the ultraviolet (in p^). As a result, as A ^ oo 
for any given solution the values for the bare mass also presented a decaying oscillatory 
behaviour 0. Clearly, these oscillating solutions cannot be obtained from the chiral limit 
procedure (i.e., mo (A) = and A — > oo) and hence cannot be chirally symmetric solutions. 



III. RESULTS 

To motivate our general arguments concerning the chiral limit, let us consider Fig. |l|, 
which shows the non-oscillating solution A{jP) and M{p'^) for supercritical coupling and 
with mo(A) = 0, for a wide range of values of A. For any given set of parameters a, mo(A), 
/i, ^, and A it is always possible to find such a non-oscillating solution. It is clear from these 
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numerical solutions that in the continuum limit (i.e., A oo) and for supercritical coupling, 
we find — > 1 for all and a mass function M{p^) which diverges proportionally to A. 
This divergent behavior of the mass was verified to the numerical accuracy of our solutions (1 
in 10^). Conversely, for subcritical coupling (i.e., a < etc), the chiral limit does exist and we 
find simply that A(/i;p2) = (p7/i2)-°5/4'r_ 

seems clear from these numerical studies that 
above critical coupling there is no finite chiral limit in the continuum quenched theory in 
any covariant gauge, even in the presence of the renormalization procedure. In the absence 
of any renormalization program, this divergent behavior of the mass is well known (see, e.g., 
Refs. 0,0) and is inevitable since in the absence of ECSB in the quenched theory there is 
only one external scale (i.e.. A). 

While the above conclusions were based on a numerical study with a specific choice of 
vertex, it is relatively straightforward to construct a general argument which applies irrespec- 
tive of this choice: Consider any vertex F which satisfies the WTI and which leads to multi- 
plicative renormalizability It automatically follows that M{p^) and y4(p^)/Z2(/i, A) are 
renormalization point independent as discussed previously. We can then define dimensionless 
quantities by appropriately scaling with A, i.e., Ji = /i/A, p'^ = p^/A^, M{p^) = M(p^)/A, 
A(p^) = y4(/i; p^)/Z2(/i, A), and niQ = mo(A)/A. Note that the renormalization condi- 
tion A{fi; fi'^) = 1 automatically determines Z2(/i,A) for a given solution for fixed A, (see 
Eq. (^)). The dimensionless functions M(p^) and A(p^) of the dimensionless variable 
can only depend on dimensionless parameters, i.e., a, ^, fl, and mg. Furthermore, since for 
any fixed A they are independent of fi (recall that we are working only in the quenched 
approximation), then it follows that they must in turn be independent of fi. Now for any 
finite A and the choice mo (A) = we have mo = 0. Hence solving for any /i and A with 
mo(A) = allows us to form M(p^) and A(p^), from which we can read off the solutions 
for A(p^) and M(p^) for any other and A with vanishing bare mass using the above rules. 
We see then that the resulting M(p^) must diverge with A as was found numerically. To 
summarize, we see that above critical coupling any vertex which satisfies the WTI and leads 
to multiplicative renormalizability will lead to a diverging mass function in the continuum 
limit for quenched QED in four dimensions. 

It is interesting to contrast this result with the well-known behavior [|I^ found in QCD, 
which leads to a well defined and finite chiral limit in the continuum. The asymptotic 
behaviour of QCD has been well-established in the ultraviolet regime due to asymptotic 
freedom in the form of a decreasing running coupling constant as{^)- This can be used in 
a renormalization group improved treatment of the ultraviolet region in the DSE study of 
DCSB in the quark propagator |[I|-p|,p!9| . 

Given that we know that the chiral limit leads to a divergent mass function above critical 
coupling, what are we then to make of the well-behaved finite solutions? Our conventional 
thinking about the chiral limit would imply that since ECSB should only increase the mass 
function above that found in its absence, then any solution above critical coupling which 
also has ECSB in the conventional sense must also correspond to a divergent mass in the 
continuum limit. Thus the finite solutions do not correspond to the chiral limit nor to any 
conventional concept of ECSB. 

It has by now probably become clear that a given set of the parameters a, mo (A), fi, ^, 
and A actually admits more than one finite solution. If the renormalization point fi is chosen 
to lie far enough in the infrared that no oscillations in the mass can occur for < /x^, then 
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specifying the renormalized mass m(/i) rather than the bare mass mo (A) leads to a unique 
solution. This corresponds to the procedure used in obtaining the numerical solutions in 
Refs. We explicitly show this behavior in Figs. |^ and ^, where for a given set of 

parameters a, mo (A), /i, ^, and A we see that there are distinct solutions. As A — oo the 
number of simultaneous solutions becomes infinite. This can be understood in the following 
way: The oscillations have a period in terms of In(p^) which is independent of A 0. Thus 
the higher is A the more oscillations can be squeezed in between the IR and the UV cut-off. 
As we increase m(/i) with everything else fixed, we push oscillations along the axis. Each 
time that a full oscillation is pushed past the UV cut-off (A^) there will be two solutions 
which give the same value for mo (A). 

It seems reasonable to expect that we can also induce these multiple solutions in QCD 
by lowering m(/i) below that which corresponds to the chiral limit. This would correspond 
to a "negative" ECSB and would then again lead to multiple solutions. A numerical test of 
this expectation is currently being pursued by implementing the renormalization program 
in a QCD-based study of the quark DSE pO|. While the DSE for the quark propagator is 
not yet well-understood in the IR, the multiple solutions should manifest themselves in the 
UV and will not depend on the detailed ansatze used for the QCD propagators and vertex 
in the IR. 



IV. SUMMARY AND CONCLUSIONS 

We have studied renormalized quenched strong-coupling QED in four dimensions in arbi- 
trary covariant gauge. We saw that there is no finite chiral limit of the renormalized theory 
on general grounds above the critical coupling. In addition, we showed that above criti- 
cal coupling for the gauge-covariance-improved treatment of the Curtis-Pennington proper 
fermion-photon vertex, there are an infinite number of renormalized finite solutions corre- 
sponding to the same bare mass in the continuum limit. All of the finite solutions have 
oscillations and differ in how far out on the momentum scale the first oscillation occurs. 

It seems likely that this behavior for the finite solutions is independent of the detailed 
vertex choice and furthermore that the same behavior can be induced in QCD by forcing the 
renormalized mass into an unphysical regime (by choosing m(/i) = M(yU^) below the value 
corresponding to the chiral limit). It also seems likely that unquenching the theory will 
not remove this rather undesirable behaviour in QED4, since the running coupling increases 
with scale rather than decreasing as in QCD. These latter conjectures are the subject of 



current investigation ||2^ . This study supports the assertion that in four dimensions strong 



coupling QED does not have a continuum limit in the conventional sense. 
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FIGURES 

2 

A(p ) for several UV cutoffs 



1.05 




M(p ) for several UV cutoffs 




FIG. 1. The behaviour of the finite renormalization A{p'^) and the mass function M{p'^) as a 
function of the ultraviolet cut-off A for mo (A) = 0. These solutions were for renormalization point 
/i^ = 10^, coupling a = 1.15, and gauge parameter ^ = 0.25. Clearly as A ^ oo we find — > 1 
for all and M{p^) diverges with A. 
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Bare vs Renormalized Mass 




FIG. 2. The relationship between the bare mass (mo(A)) and the renormalized mass (m{p)) 
for the renormalized finite solutions. These result from solving for a given m(/i) and extracting 
the corresponding mo (A). The other parameters for these solutions were renormalization point 
/x^ = 10^, coupling a = 1.25, gauge parameter ^ = 0.25 and = 10^^. The dashed horizontal line 
shows, e.g., that for mo (A) =0.1 there are three solutions. (The dashed vertical lines connecting 
the upper and lower curves are merely to guide the eye on this back-to-back log scale.) 
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Multiple solutions for a specified bare mass 
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All solutions have mj|(A) = 0.10 
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FIG. 3. The three rcnormalized finite solutions corresponding to the same bare mass, i.e., 
mo (A) = 0.1, as indicated in Fig. 2. The other parameters for these solutions were renormalization 
point //^ = 10^, coupling a = 1.25, gauge parameter ^ = 0.25 and = 10^^. (The short-long 
dashed vertical lines connecting the upper and lower curves are merely to guide the eye on this 
back-to-back log scale.) 
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